{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Lab 25 - Hypothesis testing with multiple categories\n",
    "\n",
    "We will learn how to test hypotheses when our categorical data has more than two categories.\n",
    "\n",
    "First, let's import the necessary libraries."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Jury Panels in Alameda County 2009-2010\n",
    "\n",
    "We will look at jury panel data from 2009 and 2010 collected by the American Civil Liberties Union (ACLU).  The total number of people who reported for jury duty in those years was 1,452.  See [11.2 Multiple Categories](https://www.inferentialthinking.com/chapters/11/2/Multiple_Categories.html) for more information.\n",
    "\n",
    "We can create a dataframe with this data as shown below. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create a dictionary listing each column, followed by the column values in a list\n",
    "jury_data = {\"Eligible\":[0.15, 0.18, 0.12, 0.54, 0.01],\n",
    "             \"Panels\":[0.26, 0.08, 0.08, 0.54, 0.04]}\n",
    "# pass the dictionary into the dataframe creation function as a parameter\n",
    "# also pass in labels for the rows\n",
    "jury = pd.DataFrame(data = jury_data, index = [\"Asian\",\"Black\",\"Latino/a\",\"White\",\"Other\"])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What do you think are the columns and rows of the new dataframe?  Check your answer by displaying the dataframe `jury`. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Eligible</th>\n",
       "      <th>Panels</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Asian</th>\n",
       "      <td>0.15</td>\n",
       "      <td>0.26</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Black</th>\n",
       "      <td>0.18</td>\n",
       "      <td>0.08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Latino/a</th>\n",
       "      <td>0.12</td>\n",
       "      <td>0.08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>White</th>\n",
       "      <td>0.54</td>\n",
       "      <td>0.54</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other</th>\n",
       "      <td>0.01</td>\n",
       "      <td>0.04</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          Eligible  Panels\n",
       "Asian         0.15    0.26\n",
       "Black         0.18    0.08\n",
       "Latino/a      0.12    0.08\n",
       "White         0.54    0.54\n",
       "Other         0.01    0.04"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "jury"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We will now make a bar chart of the two columns of data.  Because our dataframe `jury` already contains counts of the different categories, we do not have to use the `value_count()` function. Instead, type `jury.plot(kind = \"bar\")` below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f314ee1a898>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEaCAYAAADqqhd6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAGVNJREFUeJzt3X+cVfV95/HXW4SMBEELs24fgM40aywoYGVAjIFQtXngj4Cr1UKajezq8kgiGyVNHsFmlxC0qbXUpKhZQlIbm/graN1FxVBNk0cSFcOQgETAgpaGsdkEiUFRUSZ+9o9zZrgMM8ydmTtzZr73/Xw87mPmnPO9937mPGbec+73fM/3KCIwM7O0HFN0AWZmVnkOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBDnczcwS5HA3M0uQw93MLEHHFvXGo0aNirq6uqLe3sxsQNq4cePLEVHbWbvCwr2uro7Gxsai3t7MbECS9G/ltHO3jJlZghzuZmYJcribmSWosD53M6tuBw8epKmpiQMHDhRdSr9UU1PDmDFjGDx4cLee73A3s0I0NTVx/PHHU1dXh6Siy+lXIoK9e/fS1NREfX19t17D3TJmVogDBw4wcuRIB3s7JDFy5MgefapxuJtZYRzsHevpvnG4m5klyH3uZgmrW/xoj56/q+bDPS9i6b6ymvW01rZ23Xxxp20GDRrEhAkTWpfnzp3L4sWLmTlzJsuXL6ehoYGLLrqIe+65hxNOOKHD11myZAkzZszgggsuaL1Ac9SoUYe1Wbp0KcOGDePTn/5093+oLnC4m1nVOu6449i0adNR26xdu7bT11m2bFmlSqoYd8uYmR1FXV0dL7/8MgA33ngjp512Gu9///uZN28ey5cvB2D+/Pk88MADrc+55ZZbmDBhAlOnTmXnzp1HvOYLL7zArFmzmDx5MtOnT2f79u0Vr9vhbmZV68033+TMM89sfdx///0dtt2wYQMPPvggmzdv5rHHHjvq3FgjRoxgy5YtLFy4kOuvv/6I7QsWLOC2225j48aNLF++nE984hMV+XlKuVvGzKpWOd0yLZ588knmzJlDTU0NNTU1fOhDH+qw7bx581q/Llq06LBt+/fv56mnnuKKK65oXffWW291o/qjc7ibmVVY6TDGtkMa33nnHU444YSy/6l0l7tlzMzKcO655/Lwww9z4MAB9u/fzyOPPNJh25bunfvvv59zzjnnsG3Dhw+nvr6e1atXA9nVqJs3b654vT5yN7N+oZyhi5XW0ufeYtasWdx8883ttp0yZQqzZ89m4sSJnHTSSUyYMIERI0a02/aVV15h4sSJvOtd7+Lee+89Yvvdd9/Nxz/+cW666SYOHjzI3LlzmTRpUmV+qJwioqIvWK6GhobwzTrMeld/Hue+bds2xo0b1/PX70P79+9n2LBhvPHGG8yYMYNVq1Zx1lln9dr7tbePJG2MiIbOnusjdzOzMi1YsICtW7dy4MABrrrqql4N9p5yuJuZlemee+4puoSy+YSqmVmCygp3SbMkPS9pp6TF7WyfL2mPpE3545rKl2pmZuXqtFtG0iDgDuCPgCZgg6Q1EbG1TdP7I2JhL9RoZmZdVM6R+1RgZ0S8GBFvA/cBc3q3LDMz64lyTqiOBnaXLDcBZ7fT7nJJM4B/ARZFxO522piZtW9p+2PGu/96nU813DLlb3NzM+PGjeOuu+5i6NChFSvhG9/4Bo2Njdx+++0Ve81yVeqE6sNAXURMBB4H7mqvkaQFkholNe7Zs6dCb21m1j0tc8v87Gc/Y8iQIaxcubLokiqmnHB/CRhbsjwmX9cqIvZGRMvMN18HJrf3QhGxKiIaIqKhtra2O/WamfWK6dOnt07Pe+mllzJ58mROP/10Vq1a1dpm2LBhfO5zn2PSpElMmzaNX/7ylwDs2bOHyy+/nClTpjBlyhSefPLJI15/9erVnHHGGUyaNIkZM2b0+s9TTrhvAE6VVC9pCDAXWFPaQNLvlizOBrZVrkQzs97V3NzMY4891npXpjvvvJONGzfS2NjIihUr2Lt3LwCvv/4606ZNY/PmzcyYMYOvfe1rAFx33XUsWrSodVrga645csDgsmXLWLduHZs3b2bNmjVHbK+0TvvcI6JZ0kJgHTAIuDMinpO0DGiMiDXAJyXNBpqBXwPze7FmM7OKKJ1bZvr06Vx99dUArFixgoceegiA3bt3s2PHDkaOHMmQIUO45JJLAJg8eTKPP/44AE888QRbtx4aQPjqq6+yf//+w97r3HPPZf78+Vx55ZVcdtllvf6zlXWFakSsBda2Wbek5PsbgBsqW5qZWe9qbz7373//+zzxxBM8/fTTDB06lJkzZ3LgwAEABg8e3DqF76BBg2hubgayaXzXr19PTU1Nh++1cuVKnnnmGR599FEmT57Mxo0bGTlyZC/9ZL5C1czsMPv27ePEE09k6NChbN++nfXr13f6nA9+8IPcdtttrcvtzdX+wgsvcPbZZ7Ns2TJqa2vZvbt3BxR6bhkz6x/KGLrYF2bNmsXKlSsZN24cp512GtOmTev0OStWrODaa69l4sSJNDc3M2PGjCNG3nzmM59hx44dRATnn39+xaf4bctT/polzFP+Dmw9mfLX3TJmZglyuJuZJcjhbmaFKapbeCDo6b5xuJtZIWpqati7d68Dvh0Rwd69e486tLIzHi1jZoUYM2YMTU1NeJ6p9tXU1DBmzJhuP9/hbmaFGDx4MPX19UWXkSx3y5iZJcjhbmaWIIe7mVmCHO5mZglyuJuZJcjhbmaWIIe7mVmCHO5mZglyuJuZJcjhbmaWIIe7mVmCHO5mZglyuJuZJcjhbmaWIIe7mVmCHO5mZglyuJuZJcjhbmaWIIe7mVmCHO5mZglyuJuZJaiscJc0S9LzknZKWnyUdpdLCkkNlSvRzMy6qtNwlzQIuAO4EBgPzJM0vp12xwPXAc9UukgzM+uaco7cpwI7I+LFiHgbuA+Y0067G4G/Ag5UsD4zM+uGcsJ9NLC7ZLkpX9dK0lnA2Ih49GgvJGmBpEZJjXv27OlysWZmVp4en1CVdAxwK/BnnbWNiFUR0RARDbW1tT19azMz60A54f4SMLZkeUy+rsXxwBnA9yXtAqYBa3xS1cysOOWE+wbgVEn1koYAc4E1LRsjYl9EjIqIuoioA9YDsyOisVcqNjOzTnUa7hHRDCwE1gHbgG9HxHOSlkma3dsFmplZ1x1bTqOIWAusbbNuSQdtZ/a8LDMz6wlfoWpmliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mlqCywl3SLEnPS9opaXE72z8maYukTZJ+JGl85Us1M7NydRrukgYBdwAXAuOBee2E9z0RMSEizgRuAW6teKVmZla2co7cpwI7I+LFiHgbuA+YU9ogIl4tWXw3EJUr0czMuurYMtqMBnaXLDcBZ7dtJOla4FPAEOC8ilRnZmbdUrETqhFxR0S8B/gs8D/bayNpgaRGSY179uyp1FubmVkb5YT7S8DYkuUx+bqO3Adc2t6GiFgVEQ0R0VBbW1t+lWZm1iXlhPsG4FRJ9ZKGAHOBNaUNJJ1asngxsKNyJZqZWVd12uceEc2SFgLrgEHAnRHxnKRlQGNErAEWSroAOAi8AlzVm0WbmdnRlXNClYhYC6xts25JyffXVbguMzPrAV+hamaWIIe7mVmCHO5mZglyuJuZJcjhbmaWIIe7mVmCHO5mZglyuJuZJcjhbmaWIIe7mVmCypp+IGlLR/Tw+fsqU4eZWQX5yN3MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBJUV7pJmSXpe0k5Ji9vZ/ilJWyU9K+m7kk6pfKlmZlauTsNd0iDgDuBCYDwwT9L4Ns1+CjRExETgAeCWShdqZmblK+fIfSqwMyJejIi3gfuAOaUNIuJ7EfFGvrgeGFPZMs3MrCvKCffRwO6S5aZ8XUeuBh7rSVFmZtYzx1byxSR9BGgAPtDB9gXAAoCTTz65km9tZmYlyjlyfwkYW7I8Jl93GEkXAJ8DZkfEW+29UESsioiGiGiora3tTr1mZlaGcsJ9A3CqpHpJQ4C5wJrSBpL+APgqWbD/qvJlmplZV3Qa7hHRDCwE1gHbgG9HxHOSlkmanTf7a2AYsFrSJklrOng5MzPrA2X1uUfEWmBtm3VLSr6/oMJ1mZlZD/gKVTOzBDnczcwS5HA3M0tQRce5W3HqFj/a49fYdfPFFajEzPoDH7mbmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mlqABPf1ARS65r6lAIWZm/YyP3M3MEuRwNzNLkMPdzCxBDnczswQ53M3MEuRwNzNLkMPdzCxBDnczswQ53M3MEuRwNzNLkMPdzCxBA3puGbP29HTOoV03X1yhSsyK4yN3M7MEOdzNzBLkcDczS1BZ4S5plqTnJe2UtLid7TMk/URSs6Q/rnyZZmbWFZ2Gu6RBwB3AhcB4YJ6k8W2a/RyYD9xT6QLNzKzryhktMxXYGREvAki6D5gDbG1pEBG78m3v9EKNZmbWReV0y4wGdpcsN+XrukzSAkmNkhr37NnTnZcwM7My9OkJ1YhYFRENEdFQW1vbl29tZlZVygn3l4CxJctj8nVmZtZPlRPuG4BTJdVLGgLMBdb0bllmZtYTnYZ7RDQDC4F1wDbg2xHxnKRlkmYDSJoiqQm4AviqpOd6s2gzMzu6suaWiYi1wNo265aUfL+BrLvGzMz6AV+hamaWIIe7mVmCPOWvHbJ0RA+fv68ydRStp/sB0tkXNmD5yN3MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBDnczcwS5HA3M0uQw93MLEEOdzOzBHnKXzOzcgywqaB95G5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJcribmSXI4W5mliCHu5lZghzuZmYJKmtuGUmzgL8FBgFfj4ib22x/F/APwGRgL/AnEbGrsqWamXVf3eJHe/T8XTUVKqSPdHrkLmkQcAdwITAemCdpfJtmVwOvRMR/Ar4E/FWlCzUzs/KV0y0zFdgZES9GxNvAfcCcNm3mAHfl3z8AnC9JlSvTzMy6opxwHw3sLlluyte12yYimoF9wMhKFGhmZl3Xp/O5S1oALMgX90t6vi/fvz2CUcDL3X6BL6TzAcX7ItPj/QDeF6W8Lw6pzL44pZxG5YT7S8DYkuUx+br22jRJOhYYQXZi9TARsQpYVU5hfUVSY0Q0FF1Hf+B9kfF+OMT74pCBti/K6ZbZAJwqqV7SEGAusKZNmzXAVfn3fwz8c0RE5co0M7Ou6PTIPSKaJS0E1pENhbwzIp6TtAxojIg1wN8B35S0E/g12T8AMzMrSFl97hGxFljbZt2Sku8PAFdUtrQ+06+6iQrmfZHxfjjE++KQAbUv5N4TM7P0ePoBM7MEOdzNzBLkcDcza4ekYyS9r+g6usvhXuXyUU+ly4Mk3V1UPUWSdIakKyV9tOVRdE1FUOYjkpbkyydLmlp0XX0tIt4hm1drQOrTK1T7A0m1wH8H6ij5+SPivxVVU8HGSrohIv4yn93z28BPiy6qr0n6PDCTbHK8tWQT5f2IbLbTavMV4B3gPGAZ8BrwIDClyKIK8l1JlwP/ONCu3am60TKSngJ+CGwEftuyPiIeLKyoAuUTvN0NbAH+EFgbEV8utqq+J2kLMAn4aURMknQS8K2I+KOCS+tzkn4SEWdJ+mlE/EG+bnNETCq6tr4m6TXg3WRZ8SYgICJieKGFlaHqjtyBoRHx2aKLKJqks0oW/xb4KvAk8ANJZ0XET4qprDBvRsQ7kpolDQd+xeHTblSTg/lU3wGtn3bfKbakYkTE8UXX0F3VGO6PSLoovzCrmv1Nm+VXyLok/obsj/q8Pq+oWI2STgC+Rvapbj/wdLElFWYF8BDwHyT9BdmUIv+r2JKKkX+y/VOgPiJulDQW+N2I+HHBpXWqGrtlWj5mvQUcZAB9zLLKkzQ4Ig62WVcHDI+IZwspqh+Q9PvA+WR/H9+NiG0Fl1QISf+b/PxDRIyTdCLwTxHR788/VN1omYg4PiKOiYjjImJ4vly1wS7pi/kRa8vyiZJuKrKmPva0pP8j6WN5qBMRu6o82L8ZEdsj4o6IuD0itkn6ZtF1FeTsiLgWOAAQEa8AQ4otqTxVF+7QGmBTJc1oeRRdU4EujIjftCzkv7wXFVhPn8qncL0+X/yypA2SviTpg/nooWp0eulC3v8+uaBaijZgzz9UXbhLugb4Adksl1/Ivy4tsqaCDSoNMUnHAVUVavmR+sqIuBR4H/AwcAHwQ0k9u6vyACLphrzbcqKkVyW9li//Cvi/BZdXlLbnH34EfLHYkspTjX3uW8jG666PiDPzvsUvRsRlBZdWCEmfBT4E/H2+6r8CayLiluKqKkZ+v4L35ovPR8RBSaMjou3NaZIm6S8j4oai6+gvBur5h2oM9w0RMUXSJrL+tLckPRcRp3f65ERJupDslxfg8YhYV2Q9RZA0k+wm77vI/ojHAldFxA8KLKtPSfr9iNjeZphsqyocHgu0dkudxOEXPf68uIrKU43h/hDZ0en1ZMP9XgEGR0TV9DPbkSRtBD4cEc/ny+8F7o2IqulrlrQqIhZI+l47myMiqm14LJL+B/B54JdkFzK1jK6bWGhhZai6cC8l6QNk93v9TkS8XXQ9RZA0DbgNGEc2CmAQ8Hq1jSCS9GzbP9j21ll1ye8ud3ZEHHFP6P6uai5ikjQ8Il6V9Dslq7fkX4eR3R6wGt1OdlvE1UAD8FEO9TtXk0ZJXwe+lS//KdBYYD2FymdDrOPwrohqnGdnN7Cv6CK6o2qO3CU9EhGXSPpXsmFNKtkcEfF7BZVWqJY7upcepZbOKVIt8hFD1wLvz1f9EPhKRLxVXFXFyMe0vwfYxKH5lyIiPllcVX1L0qfyb08HTgMeJbvwEYCIuLWIurqiao7cI+KS/Gt90bX0M2/ko0Q2SboF+AVVOEQ2D/Fb80e1awDGD7RZECusZU6Zn+ePIRy6eGlA7JeqOXJvIelcYFNEvC7pI8BZwJcHwtnv3iDpFLJxzIOBRWTnIL4SETsLLayP5b8XS4FTOLwrouo+0UlaDXwyIn5RdC1Fk3RFRKzubF1/VI3h/izZ1K4TgW8AXweujIgPFFmXFUvSdrJ/bm2ngh5wJ9K6S9LDZEelxwNnAj/m8K6I2QWVVpiW6Y87W9cfVU23TInmiAhJc4DbI+LvJF1ddFF9Lb+Yq8P/7FU4SmRfRDxWdBEF+2eyT3A/IZtUr2rl135cBIyWtKJk03CguZiquqYaw/01STcA/wWYLukYqnM/XFJ0Af3M9yT9NfCPHH60Wk0X7owmm37hBuBZsvn9nwKeiohqG03272Sjpa4A/iVf10w23n1RUUV1RTV2y/xH4MPAjyPiR/mkYX8fEe8puLTCSRoF7K3GE2m+cOeQ/AR7A1nQn5M/fhMR4wstrA9JGgz8BXAN2VXLACeTTdPx522nie6Pqu6INSL+X/6H/GFJ3wL+FajG28pNA24mG99/I/BNYBRwjKSPRsR3iqyvr0XEHxZdQz9yHFn3w4j88e8cuiakWtxCdv3LKRHxGmTXygDL88d1BdZWlqo5cs8vJ5+XP14G7gc+HRGnFFpYQSQ1An9O9se7imzq3/X5JEn3Vss4d0kfiYhvlYxrPsxAGM9cKZJWkY3rfg14BlhPNsHeK4UWVgBJO4D3tv0Um88zsz0iTi2msvJV05H7drILUy5pGeYnaUD0nfWSYyPinwAkLYuI9QD5xFHFVta33p1/be9emdVx5HPIyWTTPe8AXgKagN8c9Rnpiva6JyPit5IGxO9FNYX7ZWSX2X9P0neA+zj8KtVqU3rDgTfbbBsQv7yVEBFfzb99IiKeLN2Wj32vGhExK79n6Olk/e1/Bpwh6dfA0xHx+UIL7Ftb8+7Jw6ZcyK+N2V5QTV1SNd0yLSS9G5hD1j1zHvAPwEMtR7HVQtJvgdfJ/sEdB7zRsgmoiYjBRdVWhIE8nrk3SBoDnEsW8pcAIyPihKM/Kx2SRpONnHqT7NoHyE4yHwf854Ewx3/VhXup/Ga3VwB/EhHnd9be0iPpHLIAux74Usmm4WR/xJMKKawAkj5Jti/eRzbO/amSx5aIGBC3l6skSedx6LaDWyPiu0XW0xVVHe5m+bTPM4GPAStLNr0GPBwRO4qoqwiSbiUf2+6pBwY+h7sZ2Rw7EfFvRddhVinVdELV7GjeyK9QPR2oaVlZjRcxWRqqbmpXsw7cTTYKoh74AtlViRuKLMisJ9wtY0Z2D9WImNzmpiUbImJK0bWZdYe7ZcwyLXOF/ELSxWSX3P/OUdqb9WsOd7PMTZJGkF24cxvZUMjriy3JrPvcLWPWAUnXR0TVTSpnaXC4m3VA0s8j4uSi6zDrDo+WMetYNc89ZAOcw92sY/5YawOWT6haVZP0Gu2HeMuEamYDkvvczcwS5G4ZM7MEOdzNzBLkcDczS5DD3cwsQQ53M7ME/X/KUJIaSxdBvQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "jury.plot(kind = \"bar\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "How does the distribution of eligible jurors compare with the distribution of the jury panels?\n",
    "\n",
    "To know whether the variation between the distributions is just the result of chance, we can compute a random sample from the eligible distribution and compare it with the panel distribution.\n",
    "\n",
    "First we create a variable for the eligible population and distribution."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "population = [\"Asian\",\"Black\",\"Latino/a\",\"White\",\"Other\"]\n",
    "pop_prob = [0.15, 0.18, 0.12, 0.54, 0.01]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can create a sample from this population.  The sample should be the same size as our data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['White', 'Other', 'Asian', ..., 'Asian', 'White', 'Asian'],\n",
       "      dtype='<U8')"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sample = np.random.choice(population,p = pop_prob,size = 1452)\n",
    "sample"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compute the value counts for the sample (you will have to make it a Pandas Series first).  Since we only have the probabilities of the eligible population, we want to compute the value counts as probabilities as well.  We can do this by adding the parameter `normalize = True`.  Save the probabilities of the sample in the variable `sample_probs`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "White       0.539256\n",
       "Black       0.180441\n",
       "Asian       0.146006\n",
       "Latino/a    0.121901\n",
       "Other       0.012397\n",
       "dtype: float64"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sample_probs = pd.Series(sample).value_counts(normalize = True)\n",
    "sample_probs"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<details> <summary>Answer:</summary>\n",
    "    <code>sample_probs = pd.Series(sample).value_counts(normalize = True)\n",
    "</code>\n",
    "</details>\n",
    "\n",
    "Next we will create a new column in our dataframe called `Random` that contains the probabilities from our random sample.  To do this, type `jury[\"Random\"] = sample_probs` below and run it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "jury[\"Random\"] = sample_probs"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Display the `jury` dataframe again to check that the column was added."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Eligible</th>\n",
       "      <th>Panels</th>\n",
       "      <th>Random</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Asian</th>\n",
       "      <td>0.15</td>\n",
       "      <td>0.26</td>\n",
       "      <td>0.146006</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Black</th>\n",
       "      <td>0.18</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.180441</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Latino/a</th>\n",
       "      <td>0.12</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.121901</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>White</th>\n",
       "      <td>0.54</td>\n",
       "      <td>0.54</td>\n",
       "      <td>0.539256</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other</th>\n",
       "      <td>0.01</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.012397</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          Eligible  Panels    Random\n",
       "Asian         0.15    0.26  0.146006\n",
       "Black         0.18    0.08  0.180441\n",
       "Latino/a      0.12    0.08  0.121901\n",
       "White         0.54    0.54  0.539256\n",
       "Other         0.01    0.04  0.012397"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "jury"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Plot the bar chart of the dataframe again, and the new column will be included."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f3177fece80>"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEaCAYAAADqqhd6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAHCBJREFUeJzt3X90VeWd7/H3BwQjBaFiameJmtxeq4ABgYBQC7X1x0KL4K1jB1pbtLasqkzV/ljFeq91sNNxHGs7IL2UdvxRq2itdYoV61SHrtYftASFIoqCDkqYXouoKCpK9Hv/2DvhJCbkJDnJJvt8XmudRfbez9nne84in+zz7Gc/WxGBmZnlS5+sCzAzs9JzuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7Mc2i+rFz744IOjqqoqq5c3M+uVVq9e/WJEVLbXLrNwr6qqoq6uLquXNzPrlSQ9V0w7d8uYmeWQw93MLIcc7mZmOZRZn7uZlZfdu3dTX1/Prl27si6lV6ioqGDYsGH069evU893uJtZj6ivr2fQoEFUVVUhKety9mkRwfbt26mvr6e6urpT+3C3jJn1iF27djF06FAHexEkMXTo0C59y3G4m1mPcbAXr6uflcPdzCyH3OdulmNV8+4puu3mis8Uv+MrdnSimuY6UlsxNl/1yXbb9O3bl5qamqblmTNnMm/ePE444QSuueYaamtrOe2007j11lsZMmRIm/u5/PLLmTJlCieddFLTBZkHH3xwszZXXHEFAwcO5Otf/3rn31QXONzNrMNqbqppv1Fq3ex13VhJxxxwwAGsWbNmr22WL1/e7n7mz59fqpK6jbtlzMwKVFVV8eKLLwJw5ZVXctRRR/HRj36UWbNmcc011wBwzjnn8Itf/KLpOVdffTU1NTVMmDCBTZs2vWefzzzzDFOnTmXcuHFMnjyZDRs2dPv7cLibWdl48803OfbYY5set99+e5ttV61axZ133snatWu599579zoX1uDBg1m3bh1z587l4osvfs/2OXPmsHDhQlavXs0111zDBRdcUJL3szfuljGzslFMt0yjhx56iBkzZlBRUUFFRQWnn356m21nzZrV9O8ll1zSbNvOnTt5+OGHOeuss5rWvfXWW52ovmMc7mZmXVQ4bLHlEMZ3332XIUOGFP1HpVTcLWNm1orjjz+eu+++m127drFz505+/etft9m2sXvn9ttvZ9KkSc22HXjggVRXV3PHHXcAydWna9eu7b7CUz5yN7NMFDN0sdQa+9wbTZ06lauuuqrVtuPHj2f69OmMGjWKQw45hJqaGgYPHtxq25dffplRo0ax//77s3Tp0vdsv+WWWzj//PP5zne+w+7du5k5cyajR48uzZtqgyKiW1+gLbW1teGbdZh1r+4a515TfXjRbRuHQj755JMMHz686OftC3bu3MnAgQN54403mDJlCkuWLGHs2LE99vqtfWaSVkdEbXvP9ZG7mVkb5syZwxNPPMGuXbuYPXt2jwZ7VznczczacOutt2ZdQqf5hKqZWQ4VFe6Spkp6StImSfNa2X6OpG2S1qSPL5a+VDMzK1a73TKS+gKLgJOBemCVpGUR8USLprdHxNxuqNHMzDqomCP3CcCmiHg2It4GbgNmdG9ZZmbWFcWcUD0U2FKwXA8c10q7MyVNAZ4GLomILa20MTNLXNH6mPHO76/9aYgbp/xtaGhg+PDh3HTTTQwYMKBkJdx4443U1dVx3XXXlWyfnVWqE6p3A1URMQr4LXBTa40kzZFUJ6lu27ZtJXppM7PiNM4t8/jjj9O/f38WL16cdUndpphw3wocVrA8LF3XJCK2R0TjTDg/Aca1tqOIWBIRtRFRW1lZ2Zl6zcxKYvLkyU3T855xxhmMGzeOkSNHsmTJkqY2AwcO5LLLLmP06NFMnDiRF154AYBt27Zx5plnMn78eMaPH89DDz30nv3fcccdHHPMMYwePZopU6b0zJsqUEy4rwKOlFQtqT8wE1hW2EDS3xQsTgeeLF2JZmal1dDQwL333tt0V6brr7+e1atXU1dXx4IFC9i+fTsAr7/+OhMnTmTt2rVMmTKFH//4xwBcdNFFXHLJJU3TAn/xi+8dIDh//nzuu+8+1q5dy7Jly96zvbu12+ceEQ2S5gL3AX2B6yNivaT5QF1ELAO+Imk60AC8BJzTjTWbmXVK4dwykydP5rzzzgNgwYIF3HXXXQBs2bKFjRs3MnToUPr378+0adMAGDduHL/97W8BuP/++3niiT0DBl999VV27tzZ7LWOP/54zjnnHD796U/zqU99qtvfW0tFXaEaEcuB5S3WXV7w86XApaUtzcystFqbz/13v/sd999/P4888ggDBgzghBNOYNeuXQD069evaQrfvn370tDQACTT+K5cuZKKioo2X2vx4sX88Y9/5J577mHcuHGsXr2aoUOHdtM7ey9foWpmZW3Hjh28//3vZ8CAAWzYsIGVK1e2+5xTTjmFhQsXNi23Nlf7M888w3HHHcf8+fOprKxky5aeHUDouWXMLBtFDF3sCVOnTmXx4sUMHz6co446iokTJ7b7nAULFnDhhRcyatQoGhoamDJlyntG3nzjG99g48aNRAQnnnhit0/x25Kn/DXLMU/527t1Zcpfd8uYmeWQw93MLIcc7mZmOeRwNzPLIYe7mVkOOdzNzHLI49zNLBM1N9WUdH+NQy73pnDK3+rqam6++WaGDBnS5dfevHkz06ZN4/HHH+/yvkrFR+5mVjYKp/w96KCDWLRoUdYldRuHu5mVpUmTJrF1azJ7+c6dOznxxBMZO3YsNTU1/OpXvwKSI/Lhw4fzpS99iZEjR3LKKafw5ptvArB69WpGjx7N6NGjm/2R2LVrF+eeey41NTWMGTOGFStWAMmNPM444wxOPvlkqqqquO6667j22msZM2YMEydO5KWXXirp+3O4m1nZeeedd3jggQeYPn06ABUVFdx11108+uijrFixgq997Ws0Xr2/ceNGLrzwQtavX8+QIUO48847ATj33HNZuHAha9eubbbvRYsWIYl169axdOlSZs+e3TQR2eOPP84vf/lLVq1axWWXXcaAAQN47LHHmDRpEj/96U9L+h4d7mZWNhqn/P3gBz/ICy+8wMknnwxARPCtb32LUaNGcdJJJ7F169amG3NUV1c3TRM8btw4Nm/ezCuvvMIrr7zSdBOOz33uc02v8eCDD3L22WcDcPTRR3PEEUfw9NNPA/Dxj3+cQYMGUVlZyeDBgzn99NMBqKmpYfPmzSV9rw53MysbjX3uzz33HBHR1J1yyy23sG3bNlavXs2aNWs45JBDmo62999//6bnF0772xmF++rTp0/Tcp8+fbq039Y43M2s7AwYMIAFCxbwve99j4aGBnbs2MEHPvAB+vXrx4oVK3juuef2+vwhQ4YwZMgQHnzwQSD549Bo8uTJTctPP/00zz//PEcddVT3vZk2eCikmWWimKGL3WnMmDGMGjWKpUuX8tnPfpbTTz+dmpoaamtrOfroo9t9/g033MAXvvAFJHHKKac0rb/gggs4//zzqampYb/99uPGG29sdsTeUzzlr1mOecrf3s1T/pqZWTMOdzOzHHK4m1mPyaobuDfq6mflcDezHlFRUcH27dsd8EWICLZv305FRUWn9+HRMmbWI4YNG0Z9fT3btm3LupReoaKigmHDhnX6+Q53M+sR/fr1o7q6Ousyyoa7ZczMcsjhbmaWQw53M7MccribmeWQw93MLIcc7mZmOVRUuEuaKukpSZskzdtLuzMlhaR2J7UxM7Pu0264S+oLLAJOBUYAsySNaKXdIOAi4I+lLtLMzDqmmCP3CcCmiHg2It4GbgNmtNLuSuCfgV0lrM/MzDqhmHA/FNhSsFyfrmsiaSxwWETsdfJoSXMk1Umq8yXIZmbdp8snVCX1Aa4FvtZe24hYEhG1EVFbWVnZ1Zc2M7M2FBPuW4HDCpaHpesaDQKOAX4naTMwEVjmk6pmZtkpJtxXAUdKqpbUH5gJLGvcGBE7IuLgiKiKiCpgJTA9InwPPTOzjLQb7hHRAMwF7gOeBH4eEeslzZc0vbsLNDOzjitqyt+IWA4sb7Hu8jbantD1sszMrCt8haqZWQ453M3McsjhbmaWQw53M7MccribmeWQw93MLIcc7mZmOeRwNzPLIYe7mVkOOdzNzHLI4W5mlkMOdzOzHHK4m5nlkMPdzCyHHO5mZjnkcDczyyGHu5lZDjnczcxyyOFuZpZDDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeVQUeEuaaqkpyRtkjSvle1flrRO0hpJD0oaUfpSzcysWO2Gu6S+wCLgVGAEMKuV8L41Imoi4ljgauDakldqZmZFK+bIfQKwKSKejYi3gduAGYUNIuLVgsX3AVG6Es3MrKP2K6LNocCWguV64LiWjSRdCHwV6A98oiTVmZlZp5TshGpELIqIDwHfBP53a20kzZFUJ6lu27ZtpXppMzNroZhw3wocVrA8LF3XltuAM1rbEBFLIqI2ImorKyuLr9LMzDqkmHBfBRwpqVpSf2AmsKywgaQjCxY/CWwsXYlmZtZR7fa5R0SDpLnAfUBf4PqIWC9pPlAXEcuAuZJOAnYDLwOzu7NoMzPbu2JOqBIRy4HlLdZdXvDzRSWuy8zMusBXqJqZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeWQw93MLIeKmn4gF64Y3IG2O7qvDjOzHuAjdzOzHHK4m5nlkMPdzCyHHO5mZjnkcDczyyGHu5lZDjnczcxyyOFuZpZDDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeVQUeEuaaqkpyRtkjSvle1flfSEpD9LekDSEaUv1czMitVuuEvqCywCTgVGALMkjWjR7DGgNiJGAb8Ari51oWZmVrxijtwnAJsi4tmIeBu4DZhR2CAiVkTEG+niSmBYacs0M7OOKCbcDwW2FCzXp+vach5wb1eKMjOzrtmvlDuTdDZQC3ysje1zgDkAhx9+eClf2szMChRz5L4VOKxgeVi6rhlJJwGXAdMj4q3WdhQRSyKiNiJqKysrO1OvmZkVoZhwXwUcKalaUn9gJrCssIGkMcCPSIL9r6Uv08zMOqLdcI+IBmAucB/wJPDziFgvab6k6WmzfwEGAndIWiNpWRu7MzOzHlBUn3tELAeWt1h3ecHPJ5W4LjMz6wJfoWpmlkMOdzOzHCrpUEjLn5qbaopuu272um6sxMw6wuGeE1Xz7im67earPtmNlZjZvsDdMmZmOeRwNzPLIYe7mVkOOdzNzHLI4W5mlkMOdzOzHHK4m5nlkMPdzCyHHO5mZjnkcDczy6FePf1Ahy65r+jGQszM9jE+cjczy6FefeTeXTwTopn1dj5yNzPLIYe7mVkOOdzNzHLI4W5mlkMOdzOzHHK4m5nlkMPdzCyHHO5mZjnki5jMiuSL26w3cbhb7nRozqGrPtmNlZhlx90yZmY55HA3M8shh7uZWQ4VFe6Spkp6StImSfNa2T5F0qOSGiT9benLNDOzjmg33CX1BRYBpwIjgFmSRrRo9jxwDnBrqQs0M7OOK2a0zARgU0Q8CyDpNmAG8ERjg4jYnG57txtqNDOzDiqmW+ZQYEvBcn26rsMkzZFUJ6lu27ZtndmFmZkVoUdPqEbEkoiojYjaysrKnnxpM7OyUky4bwUOK1gelq4zM7N9VDHhvgo4UlK1pP7ATGBZ95ZlZmZd0W64R0QDMBe4D3gS+HlErJc0X9J0AEnjJdUDZwE/krS+O4s2M7O9K2pumYhYDixvse7ygp9XkXTXmJnZPsBXqJqZ5ZDD3cwshzzlbzm6YnDxbasP77469gX+LCynfORuZpZDDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ453M3McsjhbmaWQw53M7MccribmeWQw93MLIcc7mZmOeRwNzPLIU/5a2bWUkemgr5iR/fV0QU+cjczyyGHu5lZDjnczcxyyOFuZpZDDnczsxxyuJuZ5ZDD3cwshxzuZmY55HA3M8shh7uZWQ55+gEzsy6ouamm6LbrZq/rxkqaKyrcJU0F/hXoC/wkIq5qsX1/4KfAOGA78HcRsbm0pZqZdV7VvHuKbru5ohsL6SHtdstI6gssAk4FRgCzJI1o0ew84OWI+J/A94F/LnWhZmZWvGL63CcAmyLi2Yh4G7gNmNGizQzgpvTnXwAnSlLpyjQzs44oJtwPBbYULNen61ptExENwA5gaCkKNDOzjuvRE6qS5gBz0sWdkp7qsddue9PBwIvNVz1e/H7P6X1fUPxZ7NFGxa18DuDPopA/iz16/LM4ophGxYT7VuCwguVh6brW2tRL2g8YTHJitZmIWAIsKaawniKpLiJqs65jX+DPIuHPYQ9/Fnv0ts+imG6ZVcCRkqol9QdmAstatFkGzE5//lvgPyMiSlemmZl1RLtH7hHRIGkucB/JUMjrI2K9pPlAXUQsA/4NuFnSJuAlkj8AZmaWkaL63CNiObC8xbrLC37eBZxV2tJ6zD7VTZQxfxYJfw57+LPYo1d9FnLviZlZ/nhuGTOzHHK4m5nlkMPdzKwVkvpI+kjWdXSWw73MpaOeCpf7Srolq3qyJOkYSZ+W9PnGR9Y1ZUGJsyVdni4fLmlC1nX1tIh4l2RerV6p7Kb8lVQJfAmoouD9R8QXsqopY4dJujQi/imd3fPnwGNZF9XTJH0bOIFkcrzlJBPlPUgy22m5+SHwLvAJYD7wGnAnMD7LojLygKQzgV/2tmt3ym60jKSHgT8Aq4F3GtdHxJ2ZFZWhdIK3W4B1wMeB5RHxg2yr6nmS1gGjgcciYrSkQ4CfRcTJGZfW4yQ9GhFjJT0WEWPSdWsjYnTWtfU0Sa8B7yPJijdJZiaIiDgw08KKUHZH7sCAiPhm1kVkTdLYgsV/BX4EPAT8XtLYiHg0m8oy82ZEvCupQdKBwF9pPu1GOdmdTvUd0PRt991sS8pGRAzKuobOKsdw/7Wk09ILs8rZ91osv0zSJfE9kl/qT/R4RdmqkzQE+DHJt7qdwCPZlpSZBcBdwAck/SPJlCL/J9uSspF+s/0sUB0RV0o6DPibiPhTxqW1qxy7ZRq/Zr0F7KYXfc2y0pPULyJ2t1hXBRwYEX/OpKh9gKSjgRNJfj8eiIgnMy4pE5L+L+n5h4gYLun9wH9ExD5//qHsRstExKCI6BMRB0TEgely2Qa7pO+mR6yNy++X9J0sa+phj0j6d0lfTkOdiNhc5sF+c0RsiIhFEXFdRDwp6eas68rIcRFxIbALICJeBvpnW1Jxyi7coSnAJkia0vjIuqYMnRoRrzQupP95T8uwnh6VTuF6cbr4A0mrJH1f0inp6KFyNLJwIe1/H5dRLVnrtecfyi7cJX0R+D3JLJf/kP57RZY1ZaxvYYhJOgAoq1BLj9QXR8QZwEeAu4GTgD9IKv6uyr2cpEvTbstRkl6V9Fq6/FfgVxmXl5WW5x8eBL6bbUnFKcc+93Uk43VXRsSxad/idyPiUxmXlglJ3wROB25IV50LLIuIq7OrKhvp/Qo+nC4+FRG7JR0aES1vTpNrkv4pIi7Nuo59RW89/1CO4b4qIsZLWkPSn/aWpPURMbLdJ+eUpFNJ/vMC/DYi7suynixIOoHkJu+bSX6JDwNmR8TvMyyrR0k6OiI2tBgm26QMh8cCTd1Sh9D8osfns6uoOOUY7neRHJ1eTDLc72WgX0SUTT+zvZek1cBnIuKpdPnDwNKIKJu+ZklLImKOpBWtbI6IKLfhsUj6e+DbwAskFzI1jq4blWlhRSi7cC8k6WMk93v9TUS8nXU9WZA0EVgIDCcZBdAXeL3cRhBJ+nPLX9jW1ll5Se8ud1xEvOee0Pu6srmISdKBEfGqpIMKVq9L/x1IcnvAcnQdyW0R7wBqgc+zp9+5nNRJ+gnws3T5s0BdhvVkKp0NsYrmXRHlOM/OFmBH1kV0RtkcuUv6dURMk/RfJMOaVLA5IuJ/ZFRaphrv6F54lFo4p0i5SEcMXQh8NF31B+CHEfFWdlVlIx3T/iFgDXvmX4qI+Ep2VfUsSV9NfxwJHAXcQ3LhIwARcW0WdXVE2Ry5R8S09N/qrGvZx7yRjhJZI+lq4C+U4RDZNMSvTR/lrhYY0dtmQSyxxjllnk8f/dlz8VKv+FzK5si9kaTjgTUR8bqks4GxwA96w9nv7iDpCJJxzP2AS0jOQfwwIjZlWlgPS/9fXAEcQfOuiLL7RifpDuArEfGXrGvJmqSzIuKO9tbti8ox3P9MMrXrKOBG4CfApyPiY1nWZdmStIHkj1vLqaB73Ym0zpJ0N8lR6SDgWOBPNO+KmJ5RaZlpnP64vXX7orLplinQEBEhaQZwXUT8m6Tzsi6qp6UXc7X5l70MR4nsiIh7sy4iY/9J8g3uUZJJ9cpWeu3HacChkhYUbDoQaMimqo4px3B/TdKlwOeAyZL6UJ6fw7SsC9jHrJD0L8AvaX60Wk4X7hxKMv3CpcCfSeb3fxh4OCLKbTTZf5OMljoLeDpd10Ay3v2SrIrqiHLslvkg8BngTxHxYDpp2A0R8aGMS8ucpIOB7eV4Is0X7uyRnmCvJQn6SenjlYgYkWlhPUhSP+AfgS+SXLUMcDjJNB3fajlN9L6o7I5YI+L/pb/In5H0M+C/gHK8rdxE4CqS8f1XAjcDBwN9JH0+In6TZX09LSI+nnUN+5ADSLofBqeP/2bPNSHl4mqS61+OiIjXILlWBrgmfVyUYW1FKZsj9/Ry8lnp40XgduDrEXFEpoVlRFId8C2SX94lJFP/rkwnSVpaLuPcJZ0dET8rGNfcTG8Yz1wqkpaQjOt+DfgjsJJkgr2XMy0sA5I2Ah9u+S02nWdmQ0QcmU1lxSunI/cNJBemTGsc5iepV/SddZP9IuI/ACTNj4iVAOnEUdlW1rPel/7b2r0yy+PIZ4/DSaZ73ghsBeqBV/b6jPyK1ronI+IdSb3i/0U5hfunSC6zXyHpN8BtNL9KtdwU3nDgzRbbesV/3lKIiB+lP94fEQ8VbkvHvpeNiJia3jN0JEl/+9eAYyS9BDwSEd/OtMCe9UTaPdlsyoX02pgNGdXUIWXTLdNI0vuAGSTdM58Afgrc1XgUWy4kvQO8TvIH7gDgjcZNQEVE9Muqtiz05vHM3UHSMOB4kpCfBgyNiCF7f1Z+SDqUZOTUmyTXPkBykvkA4H/1hjn+yy7cC6U3uz0L+LuIOLG99pY/kiaRBNjFwPcLNh1I8ks8OpPCMiDpKySfxUdIxrk/XPBYFxG94vZypSTpE+y57eATEfFAlvV0RFmHu1k67fMJwJeBxQWbXgPujoiNWdSVBUnXko5t99QDvZ/D3Yxkjp2IeC7rOsxKpZxOqJrtzRvpFaojgYrGleV4EZPlQ9lN7WrWhltIRkFUA/9AclXiqiwLMusKd8uYkdxDNSLGtbhpyaqIGJ91bWad4W4Zs0TjXCF/kfRJkkvuD9pLe7N9msPdLPEdSYNJLtxZSDIU8uJsSzLrPHfLmLVB0sURUXaTylk+ONzN2iDp+Yg4POs6zDrDo2XM2lbOcw9ZL+dwN2ubv9Zar+UTqlbWJL1G6yHeOKGaWa/kPnczsxxyt4yZWQ453M3McsjhbmaWQw53M7MccribmeXQ/wfG+aW0GxCoxAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "jury.plot(kind = \"bar\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "How does the distribution of the random sample compare to the eligible distribution?  To the panel distribution?\n",
    "\n",
    "Let's compare the panels distribution to the eligible distribution quantitatively using hypothesis testing.  We need to choose a statistic to simulate, called the *test statistic*.  In this problem, we will use something called the *Total Variation Distance (TVD)* as the test statistic.  The TVD measures the difference between two distributions.  \n",
    "\n",
    "First we will compute the TVD between the panels and eligible distribution:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Eligible</th>\n",
       "      <th>Panels</th>\n",
       "      <th>Random</th>\n",
       "      <th>Difference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Asian</th>\n",
       "      <td>0.15</td>\n",
       "      <td>0.26</td>\n",
       "      <td>0.146006</td>\n",
       "      <td>0.11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Black</th>\n",
       "      <td>0.18</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.180441</td>\n",
       "      <td>-0.10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Latino/a</th>\n",
       "      <td>0.12</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.121901</td>\n",
       "      <td>-0.04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>White</th>\n",
       "      <td>0.54</td>\n",
       "      <td>0.54</td>\n",
       "      <td>0.539256</td>\n",
       "      <td>0.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other</th>\n",
       "      <td>0.01</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.012397</td>\n",
       "      <td>0.03</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          Eligible  Panels    Random  Difference\n",
       "Asian         0.15    0.26  0.146006        0.11\n",
       "Black         0.18    0.08  0.180441       -0.10\n",
       "Latino/a      0.12    0.08  0.121901       -0.04\n",
       "White         0.54    0.54  0.539256        0.00\n",
       "Other         0.01    0.04  0.012397        0.03"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "jury[\"Difference\"] = jury[\"Panels\"] - jury[\"Eligible\"]\n",
    "jury"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What is the sum of the difference column?  \n",
    "\n",
    "To fix this, we will take the absolute differences between probabilities."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Eligible</th>\n",
       "      <th>Panels</th>\n",
       "      <th>Random</th>\n",
       "      <th>Difference</th>\n",
       "      <th>Absolute Difference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Asian</th>\n",
       "      <td>0.15</td>\n",
       "      <td>0.26</td>\n",
       "      <td>0.146006</td>\n",
       "      <td>0.11</td>\n",
       "      <td>0.11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Black</th>\n",
       "      <td>0.18</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.180441</td>\n",
       "      <td>-0.10</td>\n",
       "      <td>0.10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Latino/a</th>\n",
       "      <td>0.12</td>\n",
       "      <td>0.08</td>\n",
       "      <td>0.121901</td>\n",
       "      <td>-0.04</td>\n",
       "      <td>0.04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>White</th>\n",
       "      <td>0.54</td>\n",
       "      <td>0.54</td>\n",
       "      <td>0.539256</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Other</th>\n",
       "      <td>0.01</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.012397</td>\n",
       "      <td>0.03</td>\n",
       "      <td>0.03</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          Eligible  Panels    Random  Difference  Absolute Difference\n",
       "Asian         0.15    0.26  0.146006        0.11                 0.11\n",
       "Black         0.18    0.08  0.180441       -0.10                 0.10\n",
       "Latino/a      0.12    0.08  0.121901       -0.04                 0.04\n",
       "White         0.54    0.54  0.539256        0.00                 0.00\n",
       "Other         0.01    0.04  0.012397        0.03                 0.03"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "jury[\"Absolute Difference\"] = np.abs(jury[\"Difference\"])\n",
    "jury"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What does this do?\n",
    "\n",
    "Now take the sum of the absolute difference column.  You can use the same command as when we took the sum of a filter."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.28"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "jury[\"Absolute Difference\"].sum()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<details> <summary>Answer:</summary>\n",
    "    <code>jury[\"Absolute Difference\"].sum()\n",
    "</code>\n",
    "</details>\n",
    "\n",
    "Notice this sum is twice either the positive or negative count, so we divide it by two.  This quantity is the *total variation distance (TVD)* between the distribution of ethnicity in the eligible juror population and the panel."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.14"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_tvd = jury[\"Absolute Difference\"].sum()/2\n",
    "data_tvd"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We could have done this calculation in one line of code:\n",
    "`np.abs(jury[\"Panel\"] - jury[\"Eligible\"]).sum()/2`\n",
    "Try it below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.14"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.abs(jury[\"Panels\"] - jury[\"Eligible\"]).sum()/2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next we want to understand the distribution of the test statistic (here, the TVD) if the panels were actually from the eligible distribution.  To do this, we want to simulate random samples from the eligible distribution and compute the total variation distance between the sample and eligible distribution.\n",
    "\n",
    "First compute for one random sample, and compute the TVD between its probabilities and the eligible distribution."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0019559228650137773"
      ]
     },
     "execution_count": 74,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sample = np.random.choice(population,p = pop_prob,size = 1452)\n",
    "sample_probs = pd.Series(sample).value_counts(normalize = True)\n",
    "tvd = np.abs(sample_probs - jury[\"Eligible\"]).sum()/2\n",
    "tvd\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we want to repeat this process many times, and make a histogram of the difference TVD values.  First, use a loop to generate many samples and compute the TVD to the eligible population."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [],
   "source": [
    "tvd_list = []\n",
    "for x in range(10000):\n",
    "    sample = np.random.choice(population,p = pop_prob,size = 1452)\n",
    "    sample_probs = pd.Series(sample).value_counts(normalize = True)\n",
    "    tvd = np.abs(sample_probs - jury[\"Eligible\"]).sum()/2\n",
    "    tvd_list.append(tvd)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next make the histogram of these simulated test statistics (the TVDs)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f314e9d4940>"
      ]
     },
     "execution_count": 80,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAD9dJREFUeJzt3X2MHPV9x/H3t3YAl0ttnnKiNskZhdIaTgnyCVKlke6gJCS0wVUQIkLUSLSnqE2VqkaKk7QSqfpgKlFaKZUiq4nif9qDolIsCE2Rw6VK1ZL4eOjxIIohh4rrgEKMm6OU6qpv/7hxdLj27ePs7v32/ZJWNzP7m5nv17v3ufHs7mxkJpKkte8n+l2AJKk7DHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIdb3cmfnnntujo2NrTrmjTfe4Mwzz+xNQX00LH3C8PQ6LH3C8PQ6KH3Ozc39IDPPazSup4E+NjbGwYMHVx0zOzvL5ORkbwrqo2HpE4an12HpE4an10HpMyJeamacp1wkqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQPf2kqFoztvvBvux3Yc+1fdmvpM54hC5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhmroeekQsAD8C/hdYysyJiDgbuBsYAxaAGzLzaD1lSpIaaeUIfSoz35+ZE9X8buBAZl4EHKjmJUl90skpl+uAfdX0PmBH5+VIktrVbKAn8A8RMRcR09Wy0cw8Uk1/HxjtenWSpKZFZjYeFLE5Mw9HxLuAh4HfAvZn5qYVY45m5lknWXcamAYYHR3dPjMzs+q+FhcXGRkZaa2LNaiZPucPH+tRNW83vnljV7fnY1qeYel1UPqcmpqaW3G6+5SaCvS3rRBxO7AI/DowmZlHIuJ8YDYzL15t3YmJiTx48OCq25+dnWVycrKlmtaiZvos5UuifUzLMyy9DkqfEdFUoDc85RIRZ0bEO49PAx8GngL2AzurYTuB+9svV5LUqWbetjgK3BcRx8f/VWb+fUR8F7gnIm4FXgJuqK9MSVIjDQM9M18E3neS5a8BV9VRlCSpdX5SVJIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRDNXD5XQ6bbX6yxa3yJW5rcZre/XEMaJh6hS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVIimAz0i1kXE4xHxQDW/NSIejYhDEXF3RJxWX5mSpEZaOUL/DPDsivk7gLsy873AUeDWbhYmSWpNU4EeEVuAa4G/rOYDuBK4txqyD9hRR4GSpOZEZjYeFHEv8MfAO4HbgFuAf6mOzomIC4CHMvPSk6w7DUwDjI6Obp+ZmVl1X4uLi4yMjLTWxRrUTJ/zh4/1qJp6jW6AV95sbuz45o31FlOjYXnuwvD0Oih9Tk1NzWXmRKNxDb9TNCJ+CXg1M+ciYrLVQjJzL7AXYGJiIicnV9/E7OwsjcaUoJk+m/0ezkG3a3yJO+eb+/rahZsm6y2mRsPy3IXh6XWt9dnMb9kHgY9HxMeAM4CfAv4c2BQR6zNzCdgCHK6vTElSIw3PoWfm5zJzS2aOATcC38zMm4BHgOurYTuB+2urUpLUUCfvQ/8s8DsRcQg4B/hKd0qSJLWjuROblcycBWar6ReBy7tfkiSpHX5SVJIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKkRL13IZRmM1XZN81/hSMdc7lzQYPEKXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBWiYaBHxBkR8Z2IeDIino6IL1bLt0bEoxFxKCLujojT6i9XknQqzRyhvwVcmZnvA94PXBMRHwDuAO7KzPcCR4Fb6ytTktRIw0DPZYvV7DuqWwJXAvdWy/cBO2qpUJLUlKbOoUfEuoh4AngVeBh4AXg9M5eqIS8Dm+spUZLUjMjM5gdHbALuA34P+Fp1uoWIuAB4KDMvPck608A0wOjo6PaZmZlV97G4uMjIyEjTNdVt/vCxWrY7ugFeebOWTQ+cVnod37yx3mJqNGjP3ToNS6+D0ufU1NRcZk40Gre+lY1m5usR8Qjw88CmiFhfHaVvAQ6fYp29wF6AiYmJnJycXHUfs7OzNBrTS7fsfrCW7e4aX+LO+Zb++desVnpduGmy3mJqNGjP3ToNS69rrc9m3uVyXnVkTkRsAK4GngUeAa6vhu0E7q+rSElSY80cNp0P7IuIdSz/AbgnMx+IiGeAmYj4A+Bx4Cs11ilJaqBhoGfmvwKXnWT5i8DldRQlSWrdcJzE1ZoxVtNrFo0s7Lm2L/uVusmP/ktSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVwkCXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhWgY6BFxQUQ8EhHPRMTTEfGZavnZEfFwRDxf/Tyr/nIlSafSzBH6ErArM7cBHwB+MyK2AbuBA5l5EXCgmpck9UnDQM/MI5n5WDX9I+BZYDNwHbCvGrYP2FFXkZKkxlo6hx4RY8BlwKPAaGYeqe76PjDa1cokSS2JzGxuYMQI8C3gDzPzbyPi9czctOL+o5n5/86jR8Q0MA0wOjq6fWZmZtX9LC4uMjIy0kIL9Zo/fKyW7Y5ugFferGXTA2ct9Dq+eWPH2xi0526dhqXXQelzampqLjMnGo1rKtAj4h3AA8A3MvNPq2XPAZOZeSQizgdmM/Pi1bYzMTGRBw8eXHVfs7OzTE5ONqypV8Z2P1jLdneNL3Hn/Ppatj1o1kKvC3uu7Xgbg/bcrdOw9DoofUZEU4HezLtcAvgK8OzxMK/sB3ZW0zuB+9spVJLUHc0cNn0QuBmYj4gnqmWfB/YA90TErcBLwA31lChJakbDQM/MbwNxiruv6m45kqR2+UlRSSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFMNAlqRAGuiQVouGXREvDYGz3gx1vY9f4Ere0sZ2FPdd2vG8JPEKXpGIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFaJhoEfEVyPi1Yh4asWysyPi4Yh4vvp5Vr1lSpIaaeYI/WvANScs2w0cyMyLgAPVvCSpjxoGemb+I/DDExZfB+yrpvcBO7pclySpRZGZjQdFjAEPZOal1fzrmbmpmg7g6PH5k6w7DUwDjI6Obp+ZmVl1X4uLi4yMjLTQQr3mDx+rZbujG+CVN2vZ9MAZll7b7XN888buF1OzQfs9rcug9Dk1NTWXmRONxnV8ca7MzIg45V+FzNwL7AWYmJjIycnJVbc3OztLozG91M7Flpqxa3yJO+eH49pow9Jru30u3DTZ/WJqNmi/p3VZa322+y6XVyLifIDq56vdK0mS1I52A30/sLOa3gnc351yJEntauZti38N/DNwcUS8HBG3AnuAqyPieeAXq3lJUh81POGXmZ88xV1XdbkWSVIH/KSoJBXCQJekQhjoklSINfPm4G58ia8klcwjdEkqhIEuSYUw0CWpEAa6JBXCQJekQhjoklQIA12SCmGgS1IhDHRJKoSBLkmFWDMf/ZdK1a/LWizsubYv+1V9PEKXpEIY6JJUCANdkgphoEtSIQx0SSqEgS5JhTDQJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiG8OJc0pDq5KNiu8SVu6WB9LwxWj46O0CPimoh4LiIORcTubhUlSWpd24EeEeuAvwA+CmwDPhkR27pVmCSpNZ2ccrkcOJSZLwJExAxwHfBMNwqTVC6vAV+PTgJ9M/DvK+ZfBq7orBxJqk+rf0g6fa3guF79IYnMbG/FiOuBazLz16r5m4ErMvPTJ4ybBqar2YuB5xps+lzgB20VtbYMS58wPL0OS58wPL0OSp/vyczzGg3q5Aj9MHDBivkt1bK3ycy9wN5mNxoRBzNzooO61oRh6ROGp9dh6ROGp9e11mcn73L5LnBRRGyNiNOAG4H93SlLktSqto/QM3MpIj4NfANYB3w1M5/uWmWSpJZ09MGizPw68PUu1XJc06dn1rhh6ROGp9dh6ROGp9c11WfbL4pKkgaL13KRpEL0LNAbXSYgIk6PiLur+x+NiLEV932uWv5cRHykVzW3q91eI+KciHgkIhYj4ku9rrtVHfR5dUTMRcR89fPKXtfeqg56vTwinqhuT0bEr/S69lZ08nta3f/u6vl7W69qblcHj+lYRLy54nH9cq9rP6XMrP3G8oumLwAXAqcBTwLbThjzG8CXq+kbgbur6W3V+NOBrdV21vWi7j70eibwC8CngC/1u5ca+7wM+Olq+lLgcL/7qbHXnwTWV9PnA68enx+0Wyd9rrj/XuBvgNv63U+Nj+kY8FS/ezjZrVdH6D++TEBm/g9w/DIBK10H7Kum7wWuioiols9k5luZ+T3gULW9QdV2r5n5RmZ+G/jv3pXbtk76fDwz/6Na/jSwISJO70nV7emk1//KzKVq+RnAIL9o1cnvKRGxA/gey4/poOuo10HVq0A/2WUCNp9qTPULcAw4p8l1B0knva4l3erzE8BjmflWTXV2Q0e9RsQVEfE0MA98akXAD5q2+4yIEeCzwBd7UGc3dPr83RoRj0fEtyLiQ3UX2yyvh66+iYhLgDuAD/e7ljpl5qPAJRHxc8C+iHgoM9fC/8JacTtwV2YuDvhBbDccAd6dma9FxHbg7yLiksz8z34X1qsj9GYuE/DjMRGxHtgIvNbkuoOkk17Xko76jIgtwH3Ar2bmC7VX25muPKaZ+SywyPLrBoOokz6vAP4kIhaA3wY+X33wcFC13Wt1+vc1gMycY/lc/M/UXnETehXozVwmYD+ws5q+HvhmLr8CsR+4sXrFeStwEfCdHtXdjk56XUva7jMiNgEPArsz8596VnH7Oul1axUGRMR7gJ8FFnpTdsva7jMzP5SZY5k5BvwZ8EeZOcjv1OrkMT0vlr8Pgoi4kOVMerFHda+uh68qfwz4N5b/mn2hWvb7wMer6TNYfnX8EMuBfeGKdb9Qrfcc8NF+v5Jcc68LwA9ZPpJ7mRNeeR+kW7t9Ar8LvAE8seL2rn73U1OvN7P8IuETwGPAjn73Utdzd8U2bmfA3+XS4WP6iRMe01/udy/Hb35SVJIK4SdFJakQBrokFcJAl6RCGOiSVAgDXZIKYaBLUiEMdEkqhIEuSYX4P08ioo4aXcvzAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.Series(tvd_list).hist(density = True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Does the test statistic computed from the data look like it comes from this distribution?"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.4.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}